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Abstract 

In the present paper we focus on building simple nonperturbative analytical rela- 
tivistic models of magnetars. With this purpose in mind we first develop a method for 
generating exact interior solutions to the static and axisymmetric Einstein-Maxwell- 
hydrodynamic equations with anisotropic perfect fluid and with pure poloidal mag- 
netic field. Then using an explicit exact solution we present a simple magnetar model 
and calculate some physically interesting quantities as the surface elipticity and the 
total energy of the magnetized star. 
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1 Introduction 



The Soft Gamma Repeaters are spectacular phenomena occurring in the visible Universe. 
The giant flares detected so far show that the peak luminosities are of order 10 44 — 10 4e erg/s. 
One of the most promising and widely accepted explanations are the magnetars pQ. Mag- 
netars are believed to be neutron stars with ultra strong magnetic field responsible for the 
observed giant flares. The huge amount of energy released in the giant flares can be ex- 
plained by the existence of ultra strong magnetic fields with strength of the order (or larger 
than) 10 14 - 10 15 Gauss [21 EJ. The giant flares SGR 0526-66, SGR 1900+14 and SGR 
1806-20 detected so far reveal the existence of characteristic quasi periodic oscillations in 
the range of tenths of Hz to kHz [I] . These oscillations are believed to be seismic vibrations 
of the magnetars. If the hypotheses is true this will provide us with a tool to investigate 
the stellar interior. That is why the quasi periodic oscillations were intensively studied in 
the past years [5] [7] (and references therein). 

The study of the stelar interior by the quasi periodic oscillations require adequate models 
of the internal structure of the magnetars. In general, our understanding of magnetars as 
Soft Gamma Repeaters is intimately related to the understanding their internal structure 
and the construction of adequate models within General Relativity. Clearly, the building 
of completely realistic magnetar models is a formidable task. However, various simple 
relativistic models, more or less realistic, could be built and these models provide us with 
valuable physical insight into the internal structure of magnetars [S] - [UJ. The existing 
simple magnetar models are based on Einstein-Maxwell equations coupled to the perfect 
fluid hydro dynamical equations. In modeling magnetar equilibrium configurations two main 
approaches have been followed so far. The first approach is to numerically solve the coupled 
systems of equations [8], [9], [12], [TJ]. The second approach is perturbative - magnetar 
equilibrium configurations are studied by using perturbative techniques, i.e. the Einstein- 
Maxwell-hydrodynamic equations are solved by linearizing them about a known static and 
spherically symmetric background solution of Einstein-hydrodynamic equations and then 
expanding the perturbed equations in tensor harmonics [TP] , [TTj . [13], [15] - [T7]. Due to 
the linear character of the perturbative equations one can consider in a relatively simple 
manner more complicated magnetic field configurations as the simultaneous presence of 
poloidal and toroidal magnetic fields. 

In the present paper we also address the problem of constructing equilibrium configu- 
rations of neutron stars with ultra strong magnetic fields within the framework of General 
Relativity. Contrary to the previous approaches mentioned above, our approach here is 
fully analytical and nonperturbative and based on exact solutions. Exact solutions provide 
a route to better and deeper understanding of the inherent nonlinear character of grav- 
ity and its interaction with matter. On the other hand, the exact solutions could serve 
as tests for checking the computer codes which is important for the advent of numerical 
relativity. More precisely, in this paper we find exact interior solutions to the coupled 
Einstein-Maxwell-hydrodynamic equations describing static (nonrotating) equilibrium con- 
figurations of strongly magnetized neutron stars. The interaction of the neutron star fluid 
with the magnetic field is also taken into account to some extent. 
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2 Setting of the problem and exact solutions 



Our starting point is the coupled Einstein-Maxwell-hydrodynamic equations 

R,u = 8tt (t^ - \Tg,)j + 2 (f^FJ* - \F 2 g,)j (1) 

V V F» V = 4tt J", (2) 
V^F H = (3) 

where T^ u and Tj^ u M = [F^ a F v a — \F 2 g^ are the energy-momentum tensors of the 
neutron matter and the electromagnetic field, respectively. J M is the current which sources 
the electromagnetic field. 

Analytically solving of the coupled Einstein-Maxwell-hydrodynamic equations in the 
general case is a desperate task and therefore we need some simplifying assumptions. We 
will assume that the configurations (and the spacetime itself) are strictly static (nonrotating) 
and axially symmetric. In mathematical terms our assumptions mean that there exist one 
(hypersurface orthogonal) timelike Killing vector £ and one spacelike axial Killing vector 
rj, commuting with £ and with closed periodic orbits shrinking down to zero on the axis 
of symmetry. In adapted coordinates the Killing vectors can be written in the usual form 
£ = d/dt and £ = d/d(j) where t is the time coordinate and 4> is the azimuthal angle around 
the axis of symmetry. Our geometrical assumptions impose restrictions on the possible 
configurations of the electromagnetic field and the energy momentum tensor of the neutron 
star matter. More precisely, they require the absence of meridional convective currents and 
electric field. The geometric assumptions require also the 4-velocity of the neutron matter 
to be aligned with the timelike Killing vector £. 

The invariance of the Maxwell 2-form F under the axial Killing field rj and the absence 
of meridional currents allow us to introduce a magnetic potential $ defined by d$ = i v F. 
The Maxwell 2-form then is given by 



F = e~ 2u r] A (4) 

where e 2u = g(r], rj). The magnetic field B measured by a comoving observer with 4-velocity 
is B = i v -k F where * is the Hodge dual. 

In the models studied so far the neutron matter has been described by an isotropic 
perfect fluid with = (p + p)v^v u + pg^ u where p, p and are the energy density, the 
pressure and the 4-velocity of the fluid. The description of the neutron star matter as an 
isotropic perfect fluid is not completely satisfactory because it neglects the interaction of 
the neutron matter with the ultra strong magnetic field. This problem is highly nontrivial 
and extremely difficult to be solved completely. From first principles it is clear that the 
strong magnetic field yields anisotropy in the neutron star matter and this should be taken 
into account in the energy-momentum tensor of the matter. Indeed, since the neutron has 
anomalous magnetic momentum the neutron matter will react to the strong magnetic field 
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by polarizing itself due to the coupling of the neutron spin to the magnetic field. The things 
can get even more complicated if we take into account the possible manifestation of some 
quantum effects like the spin-spin interactions which can drive the system to some kind 
of feromagnetic-like state. In this context we should also note that the origin of the ultra 
strong magnetar magnetic fields is not completely clear and some sort of feromagnetic-like 
phase transition could give contribution. Even more, the estimated magnetic field on the 
magnetar surfaces mentioned above exceeds in fact the QED critical magnetic field value 
B c pa 10 13 G which shows that the nonlinear Euler-Heisenberg electrodynamics should be 
probably used instead of the linear Maxwell electrodynamics. The above arguments show 
that the proper description of the strong magnetic fields in the magnetars and the properties 
of the neutron matter require subtle and extremely complicated microscopic theory. The 
microscopic description of the magnetars is far beyond the scope of this paper where we are 
interested in the averaged macroscopic description which is astrophysically relevant. 

From a macroscopic point of view we can describe the interaction (response) of the 
neutron matter with (to) the ultra strong magnetic field by adding an anisotropic term to 
the energy-momentum tensor of the isotropic perfect fluid. The anisotropy will manifest 
itself in different pressures along the meridional planes of magnetic field and in transverse 
direction. Indeed, according to the statistical physics [T5], the pressure along the magnetic 
field is p = —Q while in transverse direction p tr = —Q — Bp, = p — Bp where Q is the grand 
canonical potential and p is the megnetization. In general the dependences Q(B) and p(B) 
should be highly nonlinear and can be determined only by the microscopic theory. In ultra 
strong magnetic field, as we mentioned, most of the neutron spins should be oriented in the 
direction of the magnetic field which means that p > 0. This shows that p tr < p in ultra 
strong magnetic field. As we will see later the exact solutions predict the same behaviour 
for the transverse pressure for realistic equations of state. 

The only anisotropic term which we can add and which is orthogonal to the merid- 
ional planes and consistent with the geometrical symmetries we imposed, is of the form 
ae(r))^e(r)) u where a is a scalar and e(r/) M is the unit vector along the axial Killing field rj. 
In other words we consider the following neutron star matter energy-momentum tensor 

T^p = (p + p)v^v v + pg^ u + ae(T))^e(rf)w (5) 
The energy-momentum tensor (J5j) can be also written in the form 

T = pv ® v + (p + <j)e{rj) ® e{rj) + p[g + v <g> v — e{rj) ® e{rj)) (6) 

which shows that p is the fluid pressure in the meridional planes where the magnetic field 
lays and p tr = p + a is the pressure in direction orthogonal to the meridional planes and 
therefore orthogonal to the magnetic field. 

Armed with the energy-momentum tensor (jSJ) we can write down the dimensionally 
reduced equations. Here we will perform the dimensional reduction with respect to the 
spacelike axial Killing vector r\. For this purpose we need to introduce the 3-dimensional 
Lorentzian metric 
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H = e 2u g - 7] <g> 77, 



(7) 



where e 2u = g(rj, rj). The covariant derivative associated with the metric H will be denoted 
by Di. Then for the reduced system of equations we obtain 



DiD*u = -4ire- 2u (p - p) - e Di$D % $ - Arrae 



-2u 



(8) 



n(H)ij = 8vr(p + p) Vi Vj + 8tt(p - p)e~ 2u H ij + 2D lU D 3 u + 2e- 2 "A$£ ) j$, (9) 



along with the contracted Bianchi identity 



(10) 



(p + p)D i U + Da> = -e- Au J^ + oD iU . 
Here lZ(H)ij is the Ricci tensor with respect to the 3-metric ifjj, e 



2U 



(11) 

and 

Our main task now is to solve the system of coupled partial differential equations (JH])- 
(jTTT) . Our strategy for solving (THT) (TTTT) is to "add nonlinearly" magnetic field to a known 
static and axisymmetric solution to Einstein-hydrodynamic equations (i.e. without magnetic 
field) described by the set {p°, p°, v®, u°, H^}. In order to do so we partially follow [Tj5] where 
a method for generating exact charged interior solutions was developed. We shall assume 
that u and $ depend on the space coordinates through one function x, i- e - u — u (x) an d 
$(x)- Substituting into eq.Q we find 



n(H)ij = 8tt(p + p)viVj + 8tt(p - p)e~ 2u H ij + 2 
If we impose the relations 



D iX D jX - (12) 



H, 



'j 



Hi X 



u 



p = p° e 2u ^~ 2x , 



e x-«(x) t ,o 



(13) 



and 



(14) 
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we obtain that eq.flH} is automatically satisfied since {p°,p°, v®, u°, if?-} is a solution to the 
static, axisymmetric Einstein-hydrodymanic equations by definition. Then we can use eqs. 
(ED and (dHD to find a and J^. 



a 



0^2w(x)-2 X 



du(x) \ 
dx J 



oMx) 
Aix 



dx 



(P°-P U ) 



0\_2u(x)-2x 



+ 



;4«(X) ^ 

4tt 



dx 



-2u( X ) 



dx 



AxA x . 



Ax^x, 

(15) 
(16) 



It can be checked that eq. ( ITT]) is automatically satisfied. Let us summarize the results 
in the following 

Proposition. Let {p° ,p° ,v® ,u° = X ,H^} be a solution to the Einstein-hydrodynamic 
equations with isotropic perfect fluid and u(x) is an arbitrary function of x with 

(^) < 1. Then {p^VuHij = H? j} u( X ), $(x), J<f>} 9™en by fT®, (T^, @® and W 

form a solution to the Einstein-Maxwell-hydrodynamic equations <fS))-< f77)j . 

This proposition allows to construct exact interior solutions with arbitrary equation of 
state for the background solution. The only exception is the case with stiff equation of state 
p° = p° which is very special and will not be considered here. 

The 4-dimensional metric can be easily recovered form the data we have. Namely, if 



dsf 



? x d<j) 2 + g^ 1 dx i dx j 



(17) 



is the spacetime metric of the Einstein-hydrodynamic solution, then 



ds 2 = e 2uix) d(f) 2 + e- 2u{x)+2x g^ j dx i dx j (18) 

is the spacetime metric of the Einstein-Maxwell-hydrodynamic solution, i.e of the magne- 
tized solution. 

From a physical point of view we have to impose some restrictions on the functional 
dependence u = u( X ). More precisely, in order for the new solution to possess a well defined 
axis of symmetry the function u(x) should be of the form 



u(x)=X + f(e 2x ), (19) 

where f( X ) is a regular function with /(0) = 0. In this way the new solution will inherit 
the axis of symmetry from the background solution used for its generation. 
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3 Explicit exact solution 



Now we consider a physically interesting and realistic explicit solution with a and 
vanishing on the star surface. The solution is obtained by requiring u(x) and to 
satisfy the equations of the affinely parameterized geodesies of the 2-dimensional metric 
dl 2 = du 2 + e~ 2u d§ 2 , i.e. the equations 



dx 2 



d 

dx 



-Mx) 



dx 



dx 
= 0. 



0. 



(20) 
(21) 



This requirement considerably simplifies equations (Tl5|) and ( fl6|) . The solution of the 
above equations is 



e 



2x 



(l + 6 2 e 2 x) 2 V ' 

e 2x 

= 6 tw (23) 

where 6 is an arbitrary parameteiE- One can see that u(x) is of the form (fT9l) and therefore 
the solution has a well defined axis of symmetry. The physical meaning of the parameter b 
can be uncovered as follows. For the strength of the magnetic field we have 



B 2 = -F 2 = e~ 2 * ( ^^-) g^dad.x = — -J ij d^d^. (24) 

2 WX / (1 + b 2 e 2 *f J K ' 

Taking into account the space is locally Euclidian in small neighborhood of the axis and 
the fact that e 2x \ ax i S = 0, it is not difficult to find the strength B of the magnetic field on 
the axis 

B 2 = 4b 2 . (25) 

In fact B is also the strength of the magnetic field on the north or south pole of the star 
surface. So the parameter b can be interpreted as being one half of the north pole magnetic 
field strength, i.e. b = \B Q . 

In order to be more specific we will consider a spherically symmetric background solution. 
Also we will present the background solution in the widely used Schwarzschild coordinates 
r and 9 with 



*The other parameter has been appropriately chosen in order to have a well defined axis. 
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2 2v 2-2/1 

9ee = r , g u = e x = r sin 0. 
Then our magnetized solution is as follows 



(26) 



where 



ds 2 = A 2 (gUt 2 + Qrrdr 2 + r 2 d9 2 ) + A~ 2 r 2 sin 2 9d<p 2 , (27) 

p = A- 2 p , p = A"V s (28) 

$ = ^A- 1 ^ 2 sin 2 0, (29) 

a = -~£ 2 A~ 3 (p° - p°)r 2 sin 2 9, (30) 

Jt = -B k-*{p -p°ysa?6, (31) 



A = l + l5 2 r 2 sin 2 9. (32) 



4 

The nonzero components of the magnetic field are 



B T = -B A- 1 >/gJ r cos9, (33) 
B» = B A-^. (34) 



V iJTT 

We see that when the background solution has a well defined boundary at r = R 
corresponding to the star surface where p°(R) = 0, the same is true for the magnetized 
solution, i.e. p(R) = since p = A 2 p°. Moreover, if p° also vanishes on the star surface the 
same holds for a and according to (|3"0~1) and (131]) . As we should expect the anisotropy 
pressure a is yielded by the magnetic field and vanishes for zero magnetic field. Also, as 
we discussed in Section 2, the transverse pressure p tr = p + a should be smaller than p in 
strong magnetic field. Indeed we see that for realistic equations of state for the background 
solutions, i.e. for p° > p° we have o < 0. 

In order to describe the way in which the magnetic field deforms the star we consider 
the space metric on the star surface, namely 



dl 2 = R 2 (A 2 d9 2 + AJ 2 sin 2 ddtf) , (35) 
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where A s = 1 + -^BqR sin 0. The circumference about the equator (9 — 7r/2) is 

27T cy JD 

K l m = , , (36) 



while for the polar circumference (0 = const) we have 



L p = 2 jf A SJ Rd# = ^1 + ^^i? 2 j . (37) 

The surface elipticity e sur f is given by 

£surf = z — - (38) 

Lp 

and e sur f < for B ^ 0. Therefore, for the solution under consideration the magnetic field 
elongates the star along the magnetic field - the star is prolate in shap^l. For small B 2 R 2 
we have e sur f ~ — ^BqR 2 . Here we should note that the numerical and perturbative models 
with pure poloidal magnetic field predict positive surface elipticity The reason for that 
discrepancy is the fact that the perturbative and numeric models consider the neutron star 
matter as pure isotropic perfect fluid without taking into account the anisotropy caused by 
the interaction with magnetic field. 

The next physical quantity we shall consider is the total energy M concentrated in the 
star 



M = ~ [ R\^d?x = 

47r J Star J Star 

Taking into account ( TT3"j) and ( I15p we find 



3p + a+ —B 2 

47T 



-gd 3 x. 



(39) 



M = Mo + \b 2 [ R [ 
£ J r =o Je- 

-2nA-\p - p )r 2 sin 2 0] \J \g° tt \g Q rr r 2 sin OdrdO 

where 

M = f ( Po + 3p Q ) v ^d 3 x (41) 

J Star- 
lit is worth noting that the negative surface elipticity is a characteristic of the specific solution we 
consider. In principle the more general solutions which can be generated via our proposition may give 
positive elipticity. In those cases, however, the current and anisotropy pressure a do not vanish on the 
star surface 



A" 



sm t 



+ cos 2 9 



(40) 
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is the total energy of the background solution. The explicit form of M depends of the back- 
ground solution but we can give a good approximation by using the interior Schwarzschild 
solution as a representative example of a background solution. The interior Schwarzschild 
solution is characterized by a constant energy density po = 3Mo/47ri? 3 and the metric and 
pressure are given by 



ds 2 Q = - 



2M 



1/2 



PO 



3M 
AnR 3 



'\ _ 2M „2\ V 2 



i? 3 



2M 2 
1 r 

2Mo\ 1 / 2 



1/2 



dt 2 



dr 



^ - + r 2 (^ 2 + sinW), 



(42) 



r _ (l _ ^ 



3(1 



2M \ 1 / 2 



(1 



2M ^ V2 



The mass-radius ratio satisfies the inequality 2Mq/R < 8/9. Substituting the interior 
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Schwarzschild solution into (j4"0l) and performing calculations up to terms in the order BqR 2 , 



we find 



M = M + l -BlR* (l - l -BlR^j fl-^j+O ((B R) 4 ) . (43) 

4 Conclusion 

In this paper we presented a simple method for generating exact interior solutions to the 
static and axisymmetric Einstein-Maxwell-hydrodynamic equations with anisotropic perfect 
fluid. On this base we can build simple nonperturbative analytical relativistic models of 
the magnetars. To the best of our knowledge these are the first nonperturbative analytical 
relativistic models of the magnetars with arbitrary equation of state. As an illustration 
we gave an explicit realistic exact interior solution for the magnetars and on its base we 
calculated the suraface elipticity of the star and its energy. 

The present work could be extended in several directions. It would interesting and 
important more general configurations of the magnetic field, i.e. a mixture of poloidal 
amd toroidal magnetic fields to be investigated. The next interesting extension is to add 
rotation to the star. The mentioned possible extensions are very challenging due to the 
highly nonlinear character of the Einstein equations. We hope, however, that some progress 
could be made. 
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